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ABSTRACT 

As is usual in dwarf spheroidal galaxies, today the Local Group galaxy Ursa Minor is depleted 
of its gas content. How this galaxy lost its gas is still a matter of debate. To study the history 
of gas loss in Ursa Minor, we conducted the first three-dimensional hydrodynamical simulations 
of this object, assuming that the gas loss was driven by galactic winds powered only by type II 
supernovae (SNe II). The initial gas setup and supernova (SN) rates used in our simulations are 
mainly constrained by the inferred star formation history and the observed velocity dispersion of 
Ursa Minor. After 3 Gyr of evolution, we found that the gas removal efficiency is higher when the 
SN rate is increased, and also when the initial mean gas density is lowered. The derived mass-loss 
rates are systematically higher in the central regions (< 300 pc), even though such a relationship 
has not been strictly linear in time and in terms of the galactic radius. The filamentary structures 
induced by Rayleigh-Taylor instabilities and the concentric shells related to the acoustic waves 
driven by SNe can account for the inferred mass losses from the simulations. Our results suggest 
that SNe II are able to transfer most of the gas from the central region outward to the galactic 
halo. However, other physical mechanisms must be considered in order to completely remove the 
gas at larger radii. 

Subject headings: galaxies: dwarf — galaxies: evolution — galaxies: individual(Ursa Minor) — galaxies: 
ISM — hydrodynamics — methods: numerical 


1. INTRODUCTION 


The Local Group is dynamically dominated by 
two spiral galaxies: Andromeda (M31) and the 
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Milky Way (e.g, Mateo 1998 : van den Berg 19991 


20071) . Less-luminous and less-massive galaxies re- 


side around these systems (e.s.. lMateolll998[ ). 

Some of the least luminous and least massive 
galaxies in the local universe are dwarf spheroidal 
(dSph) galaxies, both classical (Ursa Minor, Sex¬ 
tan, Draco, Fornax, Leo I, Leo H, Garina, etc) and 
ultra-faint dwarf galaxies. The classical dwarfs 
present a total mass lower than about 10® Mq and 
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an absolute V -magnitudes above —14, implying 

dynamical mass-to-lig ht ratios larger than about 


10 in solar u nits (e.g., Mateol 1998: van den Ber^ 


1999L 200^ iGre bel . Gallagher fc Harb ecld 200^ 
Grebe' 2008t Strigari et al.l2007nWalker l2013ll 
In contrast to dwarf irregular galaxies, dSphs 
are practically depleted of gas with H I col¬ 


umn d ensiti e s smaller than ^ 10^ ^ cm ^ ( Yound 
19991 2000l : Grcevich fc Putman 2009[l . domi¬ 


nated by old and intermediate-age stella r pop¬ 


ulations ( Mateol 1998 : van den BergI 2007 1. 


even 


though some recent star formation activity has 
been inferred in some dSphs (e.g., Mated 1998; 


Carrera^^^j^ 2002t Lanfranchi fc Matteuccill2004t 

DolDhin_e^_al. 2005 : Batta£liu£^_ah 200^ van den Berg 


2007: Kirby et ^ 2011: de Boer et al. 2012aibl. 

Discovered by Wilson ( 1955ll . the galaxy Ursa 
Minor is a good example of a classical dSph or¬ 
biting the Milky Way. L ocated at a heliocentric 
dista nce of about 64 kpcQII rwin fc Hatzidimitrioul 
Il995 1. Ursa Minor presents a spatial stellar distri¬ 
bution compatible with a King profile with a core 
radius of about 300 pc an d a tidal radius roughly 
between 0.9 and 1.5 kpc dlrwin fc Hatzidimitriou 


1995: Klevna et al.lll998 : Palma et al.l 20031 : Strigari et al 


20071) . Kinematic studies based on stellar ve¬ 
locity measurements show that Urs a Minor is a 


slow-rotating system (5±2 km s Mated 1998 ) 


-1 


with a velocity dispersion of about 12 km s 
inside a radius of ~36 arcmin from its nucleus 


(I Wilkinson et al.N2004l) . The integrated luminos¬ 
ity of Ursa Minor is about (2 — 3) x 10^ Lq, 
while its total mass ranges from 2 to 20 x 10^ 


about 70 and 800 ( 

Irwin & Hatzidimitriou 19951: 

Klevna et al.l 19981 

Grebel. Gallagher & Harbeck 

2003: Strigari et al. 

2007.12008 

: Wolf et al. 

2010l). 


The chemical properties of Ursa Minor are 
compatible with a relatively simple star forma¬ 
tion history characterized by a single episode of 
star formation occurring from 13 to 10 Gyr ago 
( Lanfranchi fc Matteucci l2004 l2007l) . The lack of 
recent star formation activity agrees with the in¬ 
ferred upper limits of 430 and 7000 M^t^ for the H I 
mass (lYounel [jpOOl : iGrebel. Gallagher &: HarbeckI 


20031 ). as well as the upper limit of 10° Mq for the 


^ There is a relative disagreement among the different esti¬ 
mates of the U rsa Minor’s distanc e in the literature (see the 
introduction of iPiatek et al.l2005l and references therein for 
further information). 


H II mass in Ursa Minor ( Gallagher et al.ll2003l) . 

The question that arises is what mechanism(s) 
would be r esponsible for the gas removal in 


Ursa Minor. Grebel. Gallagher fc HarbeckI ( 2003h 


listed several physical mechanisms which could 
have removed the gas component in the con¬ 
text of dSphs. Supernova (SN) feedback is one 
potential gas loss channel. SN-driven galactic 


by several authors 

(Silich & Tenorio-Tagle 19981 

Mac Low & Ferrara 

1999 

: Silich & Tenorio-Tagle| 

2001: Fragile et al. 

2003 

: Marcolini et al”^00a 

Stinson et al.l2007: 1 

Levaz et al.ll2009: Qian & Wasserburg 

2012: Ruiz et al.l 2013: Recchi] 2014). as well as in 


combination with external mechanisms such as 
ram-pressure and/or tidal stripping and ultravio¬ 


let (Uy) ba c kground radiation ( e.g., iMaver et al 
2006 . 2007t ISawala et al.l l2010f) . These works 
showed that supernovae (SNe) can heat and inject 
momentum into the interstellar medium (ISM), 
creating suitable conditions for the development 
of galactic winds. However, the efficiency of this 
process is not yet fully understood, mainly because 
of its dependence on the physical parameters of 
the galaxy (e.g., thermodynamical parameters of 
the gas, dark matter (DM) profile, etc.). 

In this work, we study the process of gas re¬ 
moval in the specific case of Ursa Minor by tak¬ 
ing into account observational estimates of its star 
formation history and by considering the effects of 
type II SNe (SNe II) feedback on the galactic mass 
loss. We aim to determine whether or not this pro¬ 
cess alone c an explain the observed properties of 
the g alaxy. Lanfranchi fc Matteuccil ( 20031 . 2004L 
20071 ). for example, were able to reproduce the 
chemical properties of this galaxy with a chem¬ 
ical evolution model characterized by a star for¬ 
mation history inferred by color-magnitude dia¬ 
grams and assuming a galactic wind triggered by 
SNe. Their galactic wind, however, is established 
in an ad-hoc manner based only on the balance 
between the thermal and the binding energies of 
the gas. Several other different mass loss processes 
have been discussed in the literature (see Section 
3.3 for further considerations on this), but we still 
lack a direct comparison to the chemical evolution 
and star formation histories of observed galaxies 
based on each of these, in particular for the case 
of Ursa Minor. 

We run three-dimensional (3D) hydrodynami- 
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cal (HD) simulations for a gas content in initial 
hydrostatic equilibrium with a static, cored DM 
gravitational field, taking into account the current 
knowledge about t he star formation history of this 
object. Following iLanfranchi &: Matteucci ( 2007 1 
and references therein, the simulations evolve the 
physical conditions of the gas during 3 Gyr, the 
interval in which stars were formed in Ursa Mi¬ 
nor. 

This paper is structured as follows. In Section 
2, we present an overview of the numerical code, as 
well as the initial conditions and numerical setup 
used in our 3D HD simulations in the context of 
the Local Group dwarf galaxy Ursa Minor. Gen¬ 
eral results on the time evolution of the gas, the 
mass-loss efficiencies in terms of the SN H rates, 
and comparisons between our results and observa¬ 
tional constraints are discussed in Section 3. The 
main conclusions obtained in this work are high¬ 
lighted in Section 4. 

2. THE HD MODEL FOR URSA MINOR 

In this section, we provide a general overview 
of the numerical code, as well as the initial con¬ 
ditions and numerical setup used in the 3D HD 
simulations of the gas content of the dSph galaxy 
Ursa Minor. 

2.1. PLUTO Code: a Brief Overview 

PLUTCd is a finite-volume / finite-difference, 
shock-capturing code designed to integrate, in one, 
two, or th ree spatial dimension s, the differential 
equations ( Mignone et al l liooi : 


dU 


= -V-T(f/)-f 5(f/), 


( 1 ) 


where U is the vector of the conservative quan¬ 
tities, T(L/') is the flux tensor, and S{U) is the 
source term. 

PLUTO can deal with the system of con¬ 
servative laws in equation o in classical HD, 
ideal/resistive magnetohydrodynamics problems, 
as well as in special relativistic hydrodynam¬ 
ics and ideal relativist i c magnetohydrodynam¬ 
ics (e.g., Mignone et al. 2007t Rossi et al. 2008t 


Burkert et al 


201T; [Bosch;_Ranionet_a]J 12012 
2012 : Mignone et al. 2012). 

In the HD case, including a radiative cooling 
function Fc and a gravitational potential $ (see 
Section 2.2 for further info rmation), U, T(U) a nd 
S{U) are defined as (e.g. iMignone et ah 200 tII 


U = 



T{U) = 


pv 

pvv -(- PI 
{EFP)v 


( 2 ) 


(3) 


and 


0 

S{U)= I -pV$ 

Fc — pv ■ V® 


(4) 


where p is the mass density, P is the thermal 
pressure, v = (ux,'Uy,Uz)^ is the fluid velocity in 
Gartesian coordinates, and I is the identity tensor 
of rank 3. 

The total energy density E is given by 


E = 


P 


p\v\ 


P- 1 


(5) 


Tesileanu et al.ll2o'^ : Porth fc Fendtll201^ Schartmann et 


considering the ideal equation of state P = (F — 
l)/oe, where F is the adiabatic index of the plasma, 
assumed as 5/3 in this work, and e = Pmt/p = 
Cs/[r(r — 1)], since P/p = Cg/F, where Pint and 
Cs are, respectively, the internal energy and the 
sound speed of the plasma. 

The Gartesian grid of the computational do¬ 
main is assumed to be fixed, i.e. invariant in time. 
Therefore, we do not take into account the ef¬ 
fects of cosmological expansion in our simulations. 
Other aspects related to the cosmological evolu¬ 
tion of the galaxy (e.g. interaction of the gas with 
cosmic background radiation, evolution of the DM 
potential, etc.) have also been neglected, as we 
consider them to be second-order processes. 

2.2. Initial Conditions and Numerical 
Setup 

-^.1. PLUTO Setup 


^http://plutocode.ph.unito.it/ 


All 3D numerical simulations performed in this 
work were made assuming a cubic domain of 3^ 
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kpc^, divided into a Cartesian grid of 256^ points, 
where the HD equations were solved using the 
supercomputer Alphacruci^ through the message 
passing interface (MPI) library for parallelization. 
The number of cores used to evolve our simulations 
in Alphacrucis varied from 256 to 520 (366 cores on 
average), which resulted in a total of about 9 x 10^ 
processor hours to run the four main simulations 
discussed in this work. 

Our HD numerical experiments evolved equa¬ 
tion m during an interval of 3 Gyr, the estimated 
duratio n of the star formation episo des in Ursa 
Minor (jhanfranchi fc Matteuccil 2004l l . 


2.2.2. Boundary Conditions 


Before studying the impact of the SNe H on the 
ISM gas of Ursa Minor, we checked the numeri¬ 
cal effects of the boundary conditions on the time 
distribution of the gas mass inside the computa¬ 
tional domain. With this purpose, we run two ad¬ 
ditional simulations using the same initial gas and 
DM configuration as in our numerical simulations 
for Mgh05SNl and Mgh05SN10 (see the next sec¬ 
tions), but with different boundary conditions: (i) 
the standard open boundary conditions, in which 
the physical quantities at the ghost cells are equal 
to those in the boundary row of the domain but 
with gradients across the boundary set as zero; 
and (ii) an outflow boundary condition, in which 
outflow velocities are constant, but inflow fluxes, 
as well as gradients across the boundaries, are set 
as zero. 


The initial gas distribution is let to evolve pas¬ 
sively for 3 Gyr, without any internal or external 
perturbations (e.g., SNe, tidal forces, etc.). As the 
gas is in hydrostatic equilibrium with DM gravi¬ 
tational potential, its total mass is supposed to be 
the same in the end of the simulation, and the ve¬ 
locity field of the gas must remain null at all times. 
However, deviations from the perfect (initial) hy¬ 
drostatic equilibrium due to th e domain discretiza¬ 
tion (e.g., IZinsrale et al.ll2nn2h can deteriorate this 
scenario, introducing spurious mass fluxes between 
consecutive mesh cells. 


In the case of the standard open boundary con¬ 
ditions, a large inflow flux of gas is observed with 
associated gas velocities reaching ~250 km s“^. 


^Cluster SGI Altix ICE 8400. Further information on the 
Alphacrucis cluster at https://lai.iag.usp.br 


These spurious fluxes produced an increment of 
1.4 X 10^^ M© inside a spherical radius of 950 pc 
(the present tidal radius of Ursa Minor). The open 
boundary condition acts as an infinite reservoir, 
which provides gas whenever the pressure equilib¬ 
rium within the domain is broken. This spurious 
inflowing flux is a major issue for the assessment of 
the time evolution of the gas mass in Ursa Minor. 

The modified outflow boundary condition be¬ 
haves in the opposite sense, with the total absence 
of a gas supply to increase the galactic mass. This 
boundary condition substantially reduced the spu¬ 
rious mass fluxes in the computational domain. 
The absolute values of the gas velocity were al¬ 
ways below of 1 km s“^. A mass of about 2.1 x 10® 
Mq was added into a spherical region of 950 pc in 
radius after 3 Gyr of evolution, which means a 
small increase of 2.7% in relation to the initial to¬ 
tal mass of gas inside the same region. This small, 
but still present, influx of gas is the result of nu¬ 
merical errors associated with domain discretiza¬ 
tion and interpolation at the boundaries. The to¬ 
tal mass inflow in this case is, however, several 
orders of magnitude smaller than that observed 
for the open boundary conditions, and small (few 
percent) compared to the initial total mass of the 
domain. Therefore, we set the boundary condi¬ 
tions for all models of Ursa Minor as outflow. 


2.2.3. The Cooling Function 

Radiative cooling for an optically thin gas is 
calculated in the P LUTO code through (e.g., 
Tesileanu et al.llioOSh 


Fr = 


dP 


= -(r-i)n2A(r), 


( 6 ) 


where n is the number density of the gas and A(T) 
is the cooling function. 

In thi s work, we adopted the cooling r ates de¬ 
rived by Wiersma. Schave fc Smith ( 2009[ l , which 
include additional effects of photo-ionization of 
heavy elements by UV and X-ray background 
radiation from galaxie s (assuming the model of 
iHaardt fc Madai]ll200lh . as well as by the cosmic 
microwave background radiation field. The cool¬ 
ing function used in our simulations was obtained 
from the in terpolation of the precomputed ta bles 
provided bv lWiersma. Schave fc Sniitb ( 2009ll . as¬ 
suming a gas number density of 0.1 cm“®, roughly 
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the initial mean number density inside the simu¬ 
lated box domain, [Fe/Hl ^ —2.13 , the median 


metallicity of Ursa Minor (jKirbv et alJl201lll , and 
a radiation field found at a redshift of about 3.8. 
Note that this redshift corresponds to an age of 
the Universe of about 1.7 Gyr, that is, roughly 
the half duration of the star formation i n Ursa 
Minor derived bv iLanfranchi fc Matteuc^ ( 2004l l 
from chemical evolution models. 

In Figure [2 we show the cooling rate as a func¬ 
tion of the temperature used in all of the numerical 
simulations performed in this work. For the sake 
of simplicity, this cooling function was used in all 
regions of the computational domain and during 
the whole computational time. 

Radiative cooling was suppressed in regions 
with temperatures lower than 1.5 x 10"^ K in 
order to mimic the equilibrium heating in the 
ISM of Ursa Minor, as well as to avoid the over 
cooling of the shocked gas produc e d by the SN 
blasts ([Mac Low fc Ferrara 19991 : Fragile et al 
I 2 OO 3 I : iFragile. Murray fc: Lin 2004 ). 


2 . 2 . 4 . SN Rates 

Even though the SN rates in galaxies are vari¬ 
able over ti me, as well as across the galaxies them¬ 
selve s le.g.. lSchmidt 19591 : Kennicuttill998l : iRecchi 
20141) . we have assumed in this work a constant 


rate in this work (in time and position) for the SN 
II explosions. 

We adopted two different time-constant SN 
rates, RsniF 1 SN per Myr and 10 SNe per Myr. 
Both rates define the range of total SNe compat- 


cn 

0) 

< 



Fig. 1.— Cooling rate function used in our nu¬ 
merical simulations of the gas loss in Ursa Minor. 


ible with those expected from chemical evolution 
models for Ursa Minor along 3 Gyr of evolution 
( Lanfranchi fc Matteuccil 20041) . In addition, the 
star formation rates related to these SN rates can 
roughly generate a total stellar mass between 0.4 
and 4x10® MgQ after 3 Gyr. This mass range is 
compatible with the inferred value of about 8x 10® 
Mrr^ in stars for U rsa Minor ( Dekel fc Woo 2003 : 
Orban et al.l 2008 ). 


2.2.5. Injection of Energy from SNe 

We assumed the following prescription to mimic 
the star formation rate, and consequently the oc¬ 
currence of SNe II in our simulations: (i) we deter¬ 
mined the SN rate timescale <RsNn = ^snii! (ii) 
for each time step, we calculated the elapsed time 
At between the occurrence of the last SN blast, 
and the time t for which calculations are 
being performed; (iii) when At > 1 computa¬ 

tional cells with number densities equal or superior 
to 0.1 cm“® are flagged as possible sites for SNe 
II take place. This is an usual ad hoc procedure 
in the literature fe.g. lKatzlll992l: ISnmmerslll993t 


Katz. W einb erg , fc He rngui^ 1996t iKav et al 

2002; Ruiz et al.n20I3l) mainly because the im¬ 
possibility of accessing molecular-cloud scales due 
to the numerical resolution of the simulations 
11.7 pc in our simulations); (iv) for each fl agged 
cell, i t is calculated the para. meter pi (e.g., iKatd 
1992 : Mori fc Umemur3l2006l) . 


Pi{t,x,y,z) = 1 - exp 




At 


ts{x,y,z) 


(7) 


which measures the probability of an i-esime star 
formation episode occurring at the time t and po¬ 
sition {x,y,z). The parameter tg is the free-fall 
timescale and C, is a dimensionless sta r formation 
rate parameter assumed to be 0.1 (e.g., Katzjll992 : 
Mori fc Umemura 2006 ). Note that the SN rates 
in our simulations are insensitive to the adopted 
value of C*, as well as to our own Equation 0, 
since i?sNii is kept fixed throughout our simula¬ 
tions. This implies that pi is only used to assist 
in the choice of the SN sites; (v) We generated 
random numbers Paieat between 0 and 1 for each 
flagged cell. If Paieat < Pi, then the associated cell 


^We have assumed a initial mass fuction of Salpeter for those 
calculations. 
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is cho sen as a possible SN site (|Mori fc Umemura 
2006ll . If there is more than one eligible site, 


then that with the highest value of pi is finally 
chosen to receive the SN ev ent, (vi) Following 
Fragile. Murray fc 113 (l2on4 . an internal energy 


of 10®^ erg is added to an approximately spherical 
volume with a two cell radius (~ 23.4 pc in our 
simulations) and centered at the elected SN site. 


where po is the initial mass density of the gas at 
the center of the galaxy, Csp is the initial sound 
speed, and x(C) = 'ycri®h(C) - 1- 

The thermal gas pressure at iteration zero is 
calculated from 


P{0 = p{0‘f- 


( 12 ) 


2.2.6. DM and Initial Gas Density Profiles 


We have assumed that the gas component 
in Ursa Minor is under the influence of grav¬ 
itational forces produced by a static, cored 
DM halo. This type of DM profile s eems to 


1995 

Salucci & Persic 1997: 

van den Bosch et al. 

2000 

de Blok & Bosma 2002 

: Klevna et al. 20031 

Simon et al.l2005l: Walker et al.l2009l: Governato et al. 

2010 

Oh et al.l201ll: Del Podo1oI 2012I: Jardel & Gebhardtl 


20121) . The DM gravitational potential <I>h 
ad opted in this work is ma thematically defined 
as ( Mac Low fc Ferrari] lOQOl ) 


$h(C)=^c 




( 8 ) 


2.2.1. The Final Setup of Our HD Models for 
Ursa Minor 

As mentioned above, the initial gas configu¬ 
ration adopted in our simulations is provided by 
equations m and m, which implies that the 
values of tq, Pc, po, and Cgg are known a priori. 
However, a precise estimate of these quantities re¬ 
quires the knowledge of the exact properties of the 
gas and DM in Ursa Minor at its very early evolu¬ 
tionary stage, which is barely known and inferred 
indirectly from photometric and spectroscopic 
data together with stellar popul ation and/or 


chemical evolution models fe.g, 

Klevna et al. 

1998. 

2003: Mated 19981 Bellazzini et al. 20021 

Garrera et al. 20021: Grebel. Gallagher & Harbeck 

2003 

: Palma et al. 

2003: iLanfranchi & Matteucci 

2004 

: Strigari et al. 

2007). 


where ^ = r/rQ and Vc^ is the maximum circular 
velocity due to this DM potential: 


= 


ro, 


(9) 


where G is the gravitational constant, and pc 
and ro are, respectively, the central mass den¬ 
sity and the characteristic radius of an isothermal, 
sph erically symmetric DM mass density profile 


Ph (Binney fc Tremainelll987 : Mac Low fc Ferrara 

199ah 


Ph(C) = 


1+e 


2 ■ 


( 10 ) 


Fo l lowing preyious works (e.g..lMa c Low fc Ferrara 
19991 : Marcolini et al. 200fil : Ruiz et aTT 2013h , we 
assumed an initial gas distribution in hydrostatic 
equilibrium with the DM grayitational potential. 
Under this condition, the initial mass density p 
for an isothermal gas follows 


P{0 = Po exp 


-r^x(0 


( 11 ) 


We have assumed ro = 300 pc and Csq = 11.5 
km“^ in all of the numerical simulations per¬ 
formed in this work. A value of 300 pc corresponds 
approximately to the core radius of the best-fitting 
King-model profile for the radial distribution of 
stars in Ursa Minor (lirwin fc Hatzidimitriouiri995t 
Klevna et al.l 1998 : Palma et al. 2003 ). A sound 
speed of 11.5 km“^ corresponds to a temperature 
of about 9544 K, which is lower than the threshold 
temperature we work out for cooling in our simu¬ 
lations (see the previous sections). This prevents 
deviations from the initial hydrostatic equilibrium 
condition. 

According to IStrigari et al.l (|2007l) . the max¬ 


imum circular velocity in Ursa Minor must be 
higher than 21 km s“^. Assuming Vc^ =21.1 km 
s“^, we used equation ([^ to estimate pc, which 
led to 6.2 X 10“^"^ g cm“^. The total mass of the 
DM halo of Ursa Minor can be obtained after in¬ 
tegrating equation (HQD in spherical coordinates, 
resulting in the following: 


Mh = in Per Q (^h - arctan^h), 


(13) 
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Table 1 

Setup model parameters of our numerical simulations for the Ursa Minor galaxy. 


Model 

-Rsnii 

(Myr-i) 

Mg 

(10^ Mo) 

Po 

(10“^^ g cm-=’) 

MghlQSNl 

1 

29.4 

4.6 

MghigSNlO 

10 

29.4 

4.6 

Mgh05SNl 

1 

7.95 

1.3 

MghOSSNlO 

10 

7.95 

1.3 


Note.— Columns are as follow: (1) model identification; (2) 
SN II rate; (3) initial mass of the gas inside the tidal radius of 
Ursa Minor; (4) initial mass density at the center of the galaxy. 
All models assume ro = 300 pc, = 3.2, = 50.0, Csq = 11.5 

km~^ (corresponding to a gas temperature of ~9544 K), Ucqq = 
21.1 km s“^, Pc — 6.2 X 10”^^ g cm~^ and — 1.51 X 10^ 
Mq (see text for further details). 





7.00 




1.75 


0.10 r 


0.05 


Fig. 2.— Spatial distribution of the gas num¬ 
ber density and pressure (from top to bottom) 
on the xy plane for the models MghlQSNl and 
MghlQSNlO at t =0 Gyr. 


Fig. 3.— Spatial distribution of the gas num¬ 
ber density and pressure (from top to bottom) 
on the xy plane for the models Mgh05SNl and 
MghOSSNlO at t =0 Gyr. 
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where = fh /rn and u is the radius of the DM 
halo. Following iMac Low fc Ferraral (199^, we 
have assumed rh = r 2 oo, the characteristic radius 
at which the mean DM density is 200 times higher 
than pcrit) the critical density of the Universe: 


?'h 


3/9c 


200pcrit 


(14) 


with pcrit = ^Hq/SttG ~ 3.7 x 10“^® g cm“^, 
where Hq is the Hubble constant, assumed to be 
equal to 71 km s“^ Mpc“^ in this work. Thus, we 
obtained = 50.0 from equation (HI, leading to 
Mh = 1.51 X 10® Mq after using equation ([13]). 

An estimate of the mass of the gas in Ursa Mi¬ 
nor at initial stages, Mg^,, is obtained from the 
integration of equation (HH over spatial coordi¬ 
nates: 


the previous models but decreasing the the value 
of Mgg/Mh from 0.1956 to 0.0526 (about a fac¬ 
tor of 4). We labeled these as Mgh05SNl and 
MghOSSNlO. A summary of the setup model pa¬ 
rameters is provided in Tabled] In Figures d] and 
|3l we present the initial spatial distribution of the 
gas number density and pressure on xy plane for 
the models listed in Tabled] 

3. RESULTS 

3.1. General Results on Time Evolution of 
the Gas Content 

In Figure [4] we show the number density, ther¬ 
mal pressure, and radial velocity maps of the gas 
on the xy central slice obtained from the model 
Mghl9SNl at the final snapshot t =3 Gyr. The 
value of the radial velocity, Urad is calculated from 





(15) 


where Ct = h/?'o and the tidal radius, rt, ^950 pc 
in the case of Ursa Minor (Irwin fc Hatzidimitriou 


199, 'ill . 


In principle, we can derive the value of po from 
equation (US if the value of Mg^ is known a pri¬ 
ori. Even though Mgp is not reliably constrained 
by observations, some reasonable assumptions can 
be made concerning its value. For example, we 
can assume that the region where the forma¬ 
tion of Ursa Minor took place in the past follows 
the primordial baryon-dark-matter ratio inferred 
from the fluctuations of the cosmic radiation back¬ 
ground temperatur e. Using the 9 yr WM AP-only 
results obtained bv iHinshaw et al.l ( 2013 1. we have 
Mgp/Mh ss 0.1956, which implies Mgp ss 2.94x 10® 
Mq, quite similar to the initial baryonic mass as¬ 
sumed in the chemical evolution models fo r Ursa 
Minor ( Lanfranchi fc Matteuccill2004 [200^ . 


Inverting equation HS) and using the values of 
Mgg, ro, F, Vc^ and Csq, we could finally obtain 
Po = 4.6 X 10“®® g cm“®. We have labeled our nu¬ 
merical simulations with po = 4.6 x 10“®® g cm“® 
and for Asnii = 1 and 10 Myr“^ as Mghl9SNl 
and Mghl9SN10, respectively. 


To check the impact of this choice on the gas 
loss driven by SNe II, we also performed two addi¬ 
tional simulations considering the same SN rates of 


■Wrad = V ■ r, (16) 

where f = r/\r\ is the unit position vector. 

The initial isothermal spherically symmetric 
distribution of the gas density and pressure is 
strongly disturbed by the SN II explosions. The 
gas content is pushed outward, spreading the ini¬ 
tial central peak seen in the density and pressure 
radial profiles. Irregular cavities are also created 
by SN blasts, as can be seen in the density plots 
(darker areas). We have confirmed that these 
low-density cavities provide low resistance and are 
natural and effic ient channels for gas flows (e.g., 
Ruiz et ^l2ni3ll . These channels are seen in the 
radial velocity map of the Figure [4] as filamentary 
structures pointing radially to the galaxy. Out¬ 
flows of gas in the central region (?;rad ^ 3 km 
s“®), as well as channels of matter flowing inward 
("^^rad ^ —3 km s“^) can be noted in the same fig¬ 
ure. 

It is important to emphasize that the re¬ 
sistance faced by the induced gas flows is re¬ 
lated to the thermal pressure of the gas sur¬ 
rounding the SN shell. SNe shells must stall 
at a radius where the ram pressure is approx¬ 
imately equal to the local thermal pressure, 
which corresponds to the condition of the sonic 
Mac h num b er of the shells reaching unity (e.g . 


Cox I972I: Ciofii 

McKee. & Bertschineed 1988; 

Ostriker & McKee 

1988 

). This happens approx- 


imately at a scale of a few tens to a hundred of 
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Fig. 4.— Spatial distribution of the gas num¬ 
ber density, thermal pressure, and radial velocity 
(from top to bottom) on the xy plane obtained 
from the model MghlQSNl at t =3 Gyr. Note 
that the scale range of the density map was re¬ 
duced by a factor of two in relation to that used 
in Figure m 


Fig. 5.— Same as Figure |4] but for model 

MghlQSNlO. 
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Fig. 6.— Same as Figure [U but for model 
MghOSSNl. Note that the scale range of the den¬ 
sity map was reduced by a factor of two in relation 
to that used in Figured] 
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Fig. 7.— Same as Figure [6] but for model 

MghOSSNlO. 
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parsec for the physical properties of the ISM in 
our model. After that, SNe remnants are expected 


to evolve as a c oustic waves through the ISM (e .g.. 
Soitzei 1982 : Falceta-Goncalves et al.l 2010al fal. 


Indeed, there is a complex distribution of shells 
with radii of a few hundreds of parsecs and sub¬ 
sonic velocities (roughly between 1 and 2 km s“^) 
seen in the radial velocity map of the Figure |4j 
indicating the presence of such acoustic waves in 
the ISM of Ursa Minor. 


Those shell-like patterns are also present in the 
pressure distribution displayed in Figure ID The 
small ring-like structure seen in Figured (mainly 
in pressure) was produced by an individual SN 
event not so far from the xy plane (more precisely, 
~ 18 pc below it). This feature has a density con¬ 
trast of about 1.5, a signature of a weak shock 
front. 


Figures El [6] and [7] are the same as Figure H] but 
referring, respectively, to the models Mghl9SN10, 
MghOSSNl, and MghOSSNlO. The features seen 
in the maps of Figure |4] are also found in those 
figures, even though they are not necessarily iden¬ 
tical. For instance, the gas content inside a radius 
of 500 pc becomes sparser and less peaked after 3 
Gyr when the SN rate is increase (e.g., compare 
the extension of the spatial distribution of the gas 
with n> 3 cm“^ in the Figures E] and El- This 
was already expected because more thermal en¬ 
ergy is injected into the gas, increasing also the 
kinetic energy available t o push it outward (e.g.. 


Mac Low fc Ferrara 19991) . The same trend is ob¬ 


served when the initial gas density is lowered (e.g., 
compare the extension of the spatial distribution 
of the gas with n > 3 cm“^ in the Figures El and 
El). This result is also expected since low-density 
environments offer less resistance to the gas mo¬ 
tions driven by the SN blasts, which is also quan¬ 
titatively corroborated by the systematic higher 
absolute values of Uj-ad shown in Figure El As it 
will be discussed in the next section, both results 
play an important role in the mass-loss process in 
Ursa Minor. 


The radial velocity maps of the gas shown in 
Figures El “El display the same features seen in the 
Urad-map of Figure m Note that the larger changes 
in the values of Urad occur in the vicinity of low- 
density cavities and/or shock fronts seen in the 
density/pressure plots. In addition, regions with 
high radial velocities (i.e., |urad| ^ 5 km s“^) are 


more abundant for the highest SN rate assumed 
in this work (i?sNii = 10 Myr“^; Figures El and El), 
as well as when the initial average gas density is 
lower (e.g., compare Figures El and El- 

In summary, we have noted two different pat¬ 
terns related to the radial velocities observed in 
our four numerical simulations: a filamentary 
structure with velocities higher than ~ 3 km s“^, 
and the superposition of concentric shells with ve¬ 
locities in the range of about 1 — 3 km s“^. A pos¬ 
sible explanation for the filamentary distribution 
of Urad is related to the Rayleigh-Taylor instabil¬ 
ity (RTI). In this process, SN explosions generate 
cavities of hot and diffuse gas which then suffer 
from the buoyancy effect. As the cavities/bubbles 
rise upward from the central region of the galaxy, 
they become stretched in filamentary morpholo¬ 
gies, similar to the ’’mushroom-like” structures 
caused by nuclear explosions. One can estimate 
the radial uprise velocity of such structures as (see 
Ruiz et aT]l2013l for details) 


VKiit) ~ O.lgt -^- , (17) 

V PISM + Pcav / 


where pisM and pcav represent the gas densities 
of the surrounding ISM and the cavities, respec¬ 
tively, and g the local acceleration of gravity from 
the DM and baryonic mass components. For the 
parameters obtained in our models, we find, at 
typical times of ~ 100 Myr, wrt "^4 — 5 km s“^, in 
absolute agreement with the velocities observed in 
the filaments. Therefore, it is probable that RTI is 
the dominant process generating these structures. 


Goncerning the concentric shells, they have ve¬ 
locities smaller than the local sound speed, i.e. at 
Mach numbers of Mg = {Sv)/cs ^ 0.1 — 0.3 (e.g., 
see Figure ED- Such behavior suggests t hat the 
shells are related to acoustic waves (e.g. , ISoitzei 


119821) , as iFalceta-Goncalves et al.l (|2010b ) had al 
ready been pointed out in the context of mas¬ 
sive galaxies in clusters. As the sound waves 
propagate radially, their amplitudes decay as does 
the wave action, resulting in a net force out¬ 
ward. The semi-analytical solution for a wave- 
driven galactic wind, assuming an isothermal and 
isotropic wind with similar parameters found in 
model Mgh05SN10, result in a specific mass-loss 
rate, M /n{r = 200 pc), equal to ^ 10“^ Mq yr“^ 
cm^ ('see lFalceta-Goncalvesll2013 ). Therefore, the¬ 
oretically, on a timescale of 1 Gyr, a total mass of 
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~ 5 X 10® Mq would have been pushed outward 
by the wave pressure in model MghOSSNlO, which 
corresponds to ~ 10% of the total mass internal to 
r = 900 pc. This value is half of the mass loss ob¬ 
served within this radius in this simulation after 1 
Gyr (see Figure|8]in Section 3.2). Combined with 
the RTI, both processes explain quantitatively well 
the mass loss inferred from this numerical simula¬ 
tion. 

3.2. Mass-loss Efficiency and SNe Rates 

In Figure m we show the time evolution of the 
total barionic mass of the gas in the ISM of Ursa 
Minor, normalized by its initial value. The relative 
masses were obtained by integrating the total gas 
mass within galactocentric radii of 0.3, 0.6, 0.95, 
and 1.5 kpc. 

For a spherical region with a radius of 1.5 kpc, 
which corresponds to almost the entire compu¬ 
tational domain, Mgh05SN10 is the only model 
presenting a considerably large mass loss (^^12%) 
after 3 Gyr. The other models present roughly 
no mass loss (model Mgh05SNl) or some accre¬ 
tion (< 3% in the worst case). We attribute the 
inferred increase of the enclosed mass in models 
Mghl9SNl and Mghl9SN10 to a combination of 
two effects. 

i. The inefficiency of blowing gas out due to the 
intrinsic higher gas densities in these models. 
Note that SN-driven perturbations advance 
to outer regions more efficiently in the case 
of the lower-density models (Mgh05SNl and 
Mgh05SN10), as shown, for instance, in the 
Figures E] “ [3 

ii. Deviations from the perfect hydrostatic 
equilibrium between gas and the DM grav¬ 
itational field due to domain discretization 
(see section 2.2). Note that spurious mass 
increments have always been smaller than 
the value (6.1%) found in our stability test 
regarding hydrostatic equilibrium discussed 
in the section 2.2.2. 

All models present gas loss for smaller radii 
(0.3, 0.6, and 0.95 kpc). It is evident that the 
higher the SN rate, the higher the mass loss. The 
same trend is observed for the gas density: lower- 
density models are prone to lose larger amounts 



Time (Gyr) 





Time (Gyr) 

Fig. 8.— From top to bottom, we show the in¬ 
stantaneous mass of the gas normalized by that 
at time zero inside a spherical radius smaller than 
1.5, 0.9, 0.6, and 0.3 kpc, respectively. Solid lines 
represent the results from model Mghl9SNl, while 
the dotted, dashed, and dashed-dotted lines refer, 
respectively, to models Mghl9SN10, Mgh05SNl, 
and MghOSSNlO. 
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galaxies fe.g. Burkert & Ruiz-Laouente 

1997 

Mac Low & Ferrarall999: Ferrara & Tolsto\ 

200 C 

Fragile. Murray & Lin 20041). However, besides 


considering a specific dSph galaxy and letting the 
SNe e xplosions d i stribu te over space and time 
(as in iRuiz et id] l201.ll l. we have also explored 


in this work how the relative mass loss is depen¬ 
dent on the radius defined to integrate it. The 
larger density of the SNe events in the the cen¬ 
tral regions of the galaxy results in the radial 
dependence of the mass loss. We also observed 
in our simulations that the gas is not completely 
removed, but pushed to for Ursa Minor’s enve¬ 
lope. Maximum mass losses of 65-85% of the 
initial total mass were obtained inside 300 pc. For 
a volume delimited by a radius of 600 pc, where 
the majority of stars in Ursa Minor has bee n de¬ 
tected (e.g., Irwin fc HatzidimitrioiJ 1995 1. the 


relative mass loss ranges from 30 to 70%. The 
gas loss decreases to about 5-50% inside 950 pc, 
the approximate tidal radius of Ursa Minor (e.g., 
Irwin fc HatzidimitripiJ 1995). 


As was already discussed in the previous sec¬ 
tion, the combinations of the RTI and the SN- 
driven acoustic waves can account for the mass 
loss inferred from our numerical simulations. 


3.3. Comparison between Numerical Sim¬ 
ulations and Observational Constraints 


Our HD simulations of the gas content in Ursa 
Minor were performed assuming a constant SN 
rate during the whole 3 Gyr. As mentioned in 
Section 2, those rates generated a total number 
of SNe at the end of the simulations compatible 
with those expecte d from chemical evolution mod - 
els for Ursa Minor ( Lanfranchi fc MatteuccilboOdI) . 
as well as a star formation history compatible with 
the inferred val ue of about 8x10 ^ Mm in stars 
for U rsa Minor ( Dekel fc Wool 2003t Orban et al 


200811 . 


Our numerical simulations showed that the 
adopted SN II rates (1 and 10 Myr“^) were able 
to transfer most of the gas from the central region 
(< 300 pc) outward to the galactic halo. However, 
SN H feedback did not completely remove the gas 
from Ursa Minor. A total mass of 1.4 x 10^ Mq re¬ 
mained inside a radius of 600 pc after 3 Gyr. This 
is two orders of magnitude higher than the upper 
limit of ^ 10® Mq for H II mass at the present time 


derived bv iGallagher et al~ ( 2003 1 from Ha obser¬ 
vations. Tentative detections of the 21 cm H I 
emission line in Ursa Minor put an upper limit of 
abou t 7000 Mm fo r the amount of neutral hydro¬ 
gen ( Young 2(300t Grebel. Gallagher fc Harbeck 


20031) . again too small compared to the remaining 


mass found from our simulations. 


This mass excess can also be noted in terms 
of the column density calculated from the spatial 
integration of the gas density along the z direction. 
After 3 Gyr, column densities inside a radius of 600 
pc from the galactic nucleus are higher than 10^^ 
cm“^ in simulations Mghl9SNl and Mghl9SN10, 


decreasing to 
^-2 


6 X 10^1 cm-2 and - 3 x lO^^ 


cm ‘ in the cases of simulations Mgh05SNl and 


Mgh05SN10, respectively. These values are higher 
than the upper limit of about 5 x 10^^ cm~^ fo r 
the H I column density derived bv I Youn j ( 2000ll . 


Therefore, our results suggest that SN H feed¬ 
back alone is insufficient to completely remove the 
ISM of Ursa Minor under the physical conditions 
adopted in this work, even though galactic winds 
have blown out a substantial amount of gas in 
our simulations (between 30 and 70% inside 600 
pc). To reconcile or at least minimize the differ¬ 
ences between our numerical results and observa¬ 
tional constraints on the gas content in Ursa Mi¬ 
nor, some additional mechanism and/or different 
initial gas/DM conditions are necessary. We list 
some possible and non-exclusive candidates in the 
next sections. 


3.3.1. Gas-to-DM Ratio 

Simulations Mgh05SNl and Mgh05SN10 started 
with ^ 8 X lO"^ Mq in gas inside 950 pc, a factor 
of ^ 3.7 lower than the initial gas mass in the 
other two simulations. It increased the mass-loss 
efficiency within 600 pc in ^23% for i?sNii = 1 
Myr“^, and ~52% for i?sNii = 10 Myr“^. These 
results suggest decreasing the gas-to-DM ratio as 
a possible way to bring the final gas mass closer 
to the observational limits for Ursa Minor. The 
main consequence of this scenario is that Ursa 
Minor could have been born in a relatively low- 
density environment in the Local Group. 

An additional possibility could be a smaller ini¬ 
tial mass of the DM halo of Ursa Minor with a 
time-dependent growth (from hierarchical accre¬ 
tion and merging of DM sub-halos) until reach- 
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ing its present-day inferred value. Such behav¬ 
ior is expected from theoretical and numerical 
works in t he context of ACDlvjl cosmologica l mod- 
els ( e.g., White &: Rees 1978 : Davis et al.l IQSi: 


Klvpin et al. 19991: Vogelsberger et al.ll2014ll . In¬ 


deed, Sawala et al.l ( 20ldl have shown (from SPH 
numerical simulations) that the time-dependent 
growth of a DM halo together with SN feed¬ 
back and an external UV radiation field can pro¬ 
duce galaxies with structural parameters similar 
to dwarf galaxies in the Local Group. 

3.3.2. SNe la 

We have ignored the feedback from SNe la in 
our simulations sin ce they are usually less frequent 
than SNe II (e .g., Lanfranchi fc Matteucci 2004 : 


Marcolini et al.l l2006[l . However, their contribu¬ 

tion to the mass-loss process may not be negligi¬ 
ble. 

SNe la occur long (1 Gyr or so) after the evo¬ 
lution of low-mass stars in binary systems, even 
during the quiescent phases of SNe II, injecting 
energy into the gas and probably contributing to 
the gas loss even after the time our simulations are 
halted. In fac^in the chemical evolution models 


by Lanfranchi fc Matteucci ( 2004 . 2007ll . the SNe 
la explosions still occur at very recent epochs, but 
at a very low rate, in Ursa Minor. They could 
also maintain Ursa Minor’s ISM at higher tem¬ 
peratures, allowing it to be less b ound to the DM 


gravitational potential w ell (e.g., iMarcolini et al 
20061 : Sawala et al.ll2010l l. 


3.3.3. Lower-density Intergalactic Medium (IGM) 

To avoid introducing a strong discontinuity to 
the gas density distribution in our simulations, we 
have used equation [Tl] beyond the tidal radius of 
Ursa Minor. It produced an external environment 
with a density ~ 10"^“® higher than typical es¬ 
timates for the IGM across the Local Group at 
present ( Mur^ 20001 : Grcevich fc Putmai] 120091 : 
Gatto et al.l I2013 I l ^These denser environments 
decelerate the galactic winds driven by the SNe 
blasts more efficiently and, as a consequence, lower 
the mass-loss efficiency. Note that this over-dense 
IGM also impacts the estimates of column densi¬ 
ties from our simulations. 


^Acronym of A Cold Dark Matter. 


3 . 3 . 4 . UV Background Radiation Field 


Some previous works have addressed the impor¬ 
tance of the UV background radiation field to the 
gas-loss process in the context of the dwarf galax- 


ies (e.g., Read. Pontzen & Vielll2006t Sawala et al. 

201c 

: Simpson et al.l2013t Milosavlievic &: Bromm 

2014 

). For example. Sawala et al. ( 

201011 and 

Simpson et al.l ( 

2013|l showed that the UV radi- 


ation field is primarily responsible for expelling 
most of the low-density intergalactic gas, even 
though the dense, cold gas in the core of the halo 
is mainly dispersed by SN feedback. 


The contribution of the UV radiation field in 
our simulations is taken into account in the radia¬ 
tive cooling function, as well as in its suppression 
below a threshold temperature (see section 2.2). 
However, the dependence on UV heating due to an 
external background has not been studied in de¬ 
tail as a free parameter in our calculations, which 
could eventually increase the mass-lo s s rates (e.g. 


Katz. We in berg . &: He rnauisd 1996 : Hoeft et al 


Falceta-Goncalvesll2ni.‘ll l . 


2006: Sawala et al.l 1201(11 : IS impson et al.l 1201,11: 


3.3.5. Tidal Stripping 


In general, tidal effects become important when 
the gravitational binding acceleration of a satel¬ 
lite galaxy is similar to the differential accelera- 
tion exerted on it due t o the host ga l axy (e .g., 
Blitz fc Robishaw 2000t Read et ^ 2006al lHl . 


This leads to the concept of the tidal radius (de¬ 
fined in section 2.2), a critical radius beyond which 
matter can be tidally stripped from the satellite 
by the host galaxy. 


iRead et ^ ( 2006al) showed through analytical 
calculations and V-body simulations that rt de¬ 
pends on the gravitational potentials of the host 
and satellite galaxies, as well as on the orbit of the 
satellite around the host and the orbit of the stars 


within the satellite . iGrebel. Gallagher & : H arbeck 
( 200.ll l. lRead et ah (2006b) and Gatto et al.n20I,'lh 


argued that tidal stripping is unimportant in¬ 
side a galactic radius of about I kpc for most 
of the Local Group dSph galaxies observed up 
to now. However, these calculations are based 
only on dynamical considerations, ruling out 
any non-gravitational energy release inside the 
tidal radius influence. For example, the en¬ 
ergy released by SNe can weaken the gravita- 
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tional potential, changing the DM density pro¬ 
files from cu sp to core during the galaxy’s evo- 


lution fe.g. 

Mashchenko. Couchman & Wadslev 

2006 

: Pasetto et al.l 2010l 

Pontzen & Governato 


2013. 


The effects of the tidal stripping mechanism on 
the gas removal in Ursa Minor will be investigated 
in future works. 

3.3.6. Ram-pressure Stripping 


Galaxies moving through the IG M are subject 


to the ram-pressure stripping if ( 

Gunn & Gott 

1972: Grebel. Gallagher & Harbeck 

200.31): 


PlGM^gal ^ Pgal^gab 


(18) 


where pigm is the intergalactic gas density, Ugai 
is the galactic velocity dispersion, and pgai and 
fgai are, respectively, the mean density and the 
velocity of the galaxy in relation to IGM. This 
simple analytical expression i s supported by more 


detai l ed HD simulations le.g . iMcGonnachie et ah 
20071 : iMaver et al.1 200fil: see iGatto et al.l 201.0l for 
some criticism), even though it is formally valid 


only for a d imensionless galaxy fsee IPasetto et al 


2 OI 2 II 2 OI 5 I for a generalization of equation 
the case of a galaxy with a non-zero size). 

In the case of Ursa Minor, the mean gas num¬ 
ber density (mean column density inside 600 
pc) must be smaller than about 2 x 10“^ cm“^ 
(7 X 10^® cm“®) in order for ram-pressure strip¬ 
ping to remove efficiently its gas component. We 
have adopted in these conservative calculations 
PiGM « 1.3 X 10“^ cm“^, corresponding to the 
minimum average particle density expected at dis- 
tances smaller than 90 kpc f rom the Mil k y Wa y 
( Grcevich fc PutmanI 2000 : Gatto et al. 20131) . 


"Cgal 

CTgal 


(Piatek et al. 


a 162 km s 
12 km s“^ (Wilkinson et al 


.Jpl), 

2004 . 


and 


We can realize that the derived upper limit of 
2 X 10“^ cm“^ is at least a factor of 26 smaller 
than the mean number density found in our nu¬ 
merical simulations after 3 Gyr of evolution. At a 
first glance, this indicates that ram-pressure strip¬ 
ping could play a small role in removing gas from 
Ursa Minor. However, Equation (ITSl) does not take 
into account nonlinear effects behind gas-shocking 
interactions, which could cha nge our naive esti - 
mate for that critical density ( Gatto et al 
Besides, a decrease in the initial gas-to-DM ratio 


could increase the amount of gas stripped by ram- 
pressure effects. 

We will address the effect of the ram-pressure 
stripping mechanism in the context of Ursa Minor 
in a forthcoming work. 


4. CONCLUSIONS 


In this work, we presented the results from 3D 
HD simulations of the gas content of the dwarf 
galaxy Ursa Minor, emulating 3 Gyr of its evo¬ 
lution. So far, this is the first direct attempt to 
model the gas mass evolution of a particular dSph 
galaxy of the Local Group in terms of 3D HD sim¬ 
ulations. We used a computational grid of 256^ 
points, distributed uniformly across a box of 3^ 
kpc^, resulting in a spatial resolution of about 11.7 
pc. 


Initially, the isothermal gas is placed in hydro¬ 
static equilibrium with a cored, static DM gravi¬ 
tational field. The total DM mass was constrained 
by the minimum value of the total velocity disper¬ 
sion of Ursa Minor. We also assumed two differ¬ 
ent values for the gas-to-DM mass ratio: about 
0.20 and 0.05, the former compatible with th e 
9 yr WMAP-only results ( Hinshaw et al.l 2013 ). 
These initial equilibrium configurations were dis¬ 
turbed by SN H explosions, at constant rates of 1 
and 10 Myr“^, in agreement with the estimates 
from chemical evolution models for this object 
( Lanfranchi fc Matteucci 2004 2007). 


From the star formation history for the dSph 
galaxy Ursa Minor assumed in this work, we con¬ 
clude the following. 


i. Gas is spread outward, erasing the initial 
central peak in the density and pressure dis¬ 
tributions. Irregular, low-density cavities 
were also created by SN-driven shocks. The 
highest spatially averaged radial velocities 
were found in the case of the highest SN rate 
assumed in this work. 

ii. Filamentary structures seen in the radial ve¬ 
locity maps (urad ^ 3 km s“^) and possibly 
induced by the RTI are predominantly asso¬ 
ciated with gas outflows. Some channels of 
matter flowing inward (urad ~3 km s“^) 
are also seen in our simulations, as well as 
a complex distribution of shell-like patterns 
with subsonic velocities roughly between 1 
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and 2 km s ^, associated with the propaga- in the context of Ursa Minor galaxy in a forthcom- 

tion of acoustic waves through the ISM. ing work. 


iii. The efficiency of the gas removal is higher 
when the SN rate is increased, as expected 
because more thermal energy is injected into 
the gas. The same trend is observed when 
the initial gas density is lowered, since low- 
density environments offer less resistance to 
the gas motions and the energy is redis¬ 
tributed over less mass, resulting in larger 
outward velocities. 


iv. The induced gas loss by SNe is differential in 
terms of the galactic radius and variable in 
time, independent of the model parameters 
considered in this work (see Figure |S]). After 
3 Gyr, the derived mass-loss rates reached 
their maximum inside 300 pc, with values 
ranging from about 65 to 85% of the initial 
mass depending on the assumed initial gas 
density and SN II rate. For a spherical ra¬ 
dius smaller than 600 pc, where the major¬ 
ity of stars in Ursa Minor has been detected 
(e.g., Irwin fc HatzidimitrioulHo^ l. the rel¬ 
ative mass-loss rates ranged from 30 to 70%. 
The gas loss decreases to about 5%-50% in¬ 
side 950 pc, roughly the tidal radius of Ursa 
Minor. 
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